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Abstract 

Using an effective theory approach, we calculate the neutron electric dipole moment (nEDM) 
in the minimal left-right symmetric model with both explicit and spontaneous CP violations. We 
integrate out heavy particles to obtain flavor-neutral CP-violating effective Lagrangian. We run the 
Wilson coefficients from the electroweak scale to the hadronic scale using one-loop renormalization 
group equations. Using the state-of-the-art hadronic matrix elements, we obtain the nEDM as a 
function of right-handed W-boson mass and CP-violating parameters. We use the current limit on 
nEDM combined with the kaon-decay parameter e to provide the most stringent constraint yet on 
the left-right symmetric scale M\y R > (10 db 3) TeV. 
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I. INTRODUCTION 



The electric dipole moment (EDM) of an elementary particle has been a subject of strong 
interest from both experimentalists and theorists for over half a century [H, 0]. A non- 
vanishing intrinsic EDM indicates violation of time- reversal (T), parity (P) and charge- 
conjugation-parity (CP) invariance at fundamental level. Although the standard model (SM) 
of particle physics predicts non-vanishing EDM for leptons and quarks from both the CP 
phase in the Cabbibo-Kobayashi-Moskawa (CKM) matrix and the quantum chromodynamic 
(QCD) #-angle, it is widely believed that new physics must exit, allowing EDMs at least 
competitive with or even dominate the SM predictions 0, 0]. New CP-violating physics 
is needed to explain, for example, the baryon number asymmetry observed in the universe 
today 0]. In this article, we are interested in the EDM of a strongly-interacting bound 
state — the free neutron. Because of its charge-neutrality, the neutron EDM (nEDM) is 
relatively "easier" to measure. The current bound is 2.9 x 10~ 26 ecm j5|, which is already 
very constraining for new physics. The upcoming experiments will enhance the current 
sensitivity by one to two orders of magnitude 0], which either will rule out many of the 
new physics models under investigation, or will provide the first opportunity to reveal an 
intrinsic nEDM. 

Critical to understanding the experimental data is a theoretical nEDM calculation with 
controlled precision of non-perturbative strong-interaction physics. The result is indispens- 
able for extracting or constraining the new interaction parameters. In the past, many 
calculations have been made in the literature, and most of which are done in the context 
of ad hoc models 0, H, @, 0, 11, 12, 13]. In one class of models, the quark EDM is first 



obtained, and the neutron EDM is calculated through constituent quark models. In an- 
other class of models, the T-odd pion-nucleon interaction vertices are first derived and then 
hadron physics effect is calculated through pion loops. The relationship among different 
contributions is often unclear and confusing. Depending on different modelings, there are 
often large uncertainties in the final result. 

In this paper, we follow an effective theory approach to calculate the nEDM in the left- 
right symmetric model (LRSM) [141 ] . The model was motivated by the hypothesis that 
parity is a perfect symmetry at high-energy, and is broken spontaneously at low-energy due 
to the asymmetric vacuum. This model has a number of attractive features, including a 
natural explanation of weak hyper-change in terms of baryon and lepton numbers, existence 
of right-handed neutrinos and the seesaw mechanism for neutrino masses, and possibility 
of spontaneous CP violation. In a recent paper, we found a complete solution of the CP 
violation structure of the minimal left-right symmetric model (mLRSM) [15(. Our goal 
here is to derive a factorization formula for nEDM in this model, with QCD and other 
short- distance physics in the Wilson coefficients, and with long-distance physics in hadronic 
matrix elements ready for, for example, lattice QCD calculations. Using the state-of-the- 
art hadronic matrix elements, we derive the best constraints on the model parameters. In 
particular, we find the most stringent bound yet on the left-right symmetric scale 10 ± 3 
TeV, which is beyond the detection capability of the Large Hadron Collider (LHC) fl6| . 

Before starting, let us make a number of relevant comments. First of all, it is possible that 
the entire nEDM to be measured can be explained by the so-called QCD #-term, a term in the 
QCD Lagrangian which will contribute to nEDM due to instanton effect. Its contribution to 
nEDM has been calculated in several different ways [TtJ ■ The nEDM constraint on #-term is 
so strong 9 < 10~ 10 [17] that there should be some mechanism, for example the Peccei-Quinn 
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symmetry [18], to cure this so-called strong CP problem. We will not consider the induced 
^-contribution to nEDM in mLRSM. However, this term seems unlikely to be the only or 
the most important source for nEDM. A typical beyond-SM physics model allow natural 
sizes of the nEDM on the order currently been probed by experiments. Second, there has 
been a number of papers in the literature about the perturbative QCD effects for nEDM in 



various versions of LRSM without general CP structure |10l . 1 191 . |20|. We will use some of 
these results to make a coherent formulation, taking into account various effects consistently. 
Finally, the EDM for spin s has an interaction term in the hamiltonian H = —ds- E/\s\, 
which corresponds to the following term in the effective lagrangian density 

C = -~d E ^ v i l5 if)F>" , (1) 

where a^ u = ^[j^jju] and is the electromagnetic field strength and if) is the spin-1/2 
Dirac field. 

The presentation of the paper is organized as follows: In Sec. II, we collect all the P-odd 
and CP-odd operators up to dimension-six and discuss their short- distance QCD effects. In 
Sec. Ill, we calculate the Wilson coefficients of quark EDM, CDM operators, four-quark 
operators and the Weinberg operator in the framework of mLRSM using effective theory 
by integrating out heavy particles. In Sec. IV, we study nEDM in mLRSM numerically 
and get the lower bound of the mass of the righthanded VT-boson (Wr). It turns out that 
the biggest contribution comes from four-quark operators. This property is true not only 
in mLRSM, but for a large category of models, for example two-Higgs-doublet Models (see 



Ref. [21] for a good review). Recently, we have made a comprehensive study of the matrix 



elements of the CP-odd and P-odd four-quark operators in the neutron state 22j, which 



makes it possible to calculate the nEDM more reliably. We conclude the paper in Sec. V. 



II. GENERAL CP- VIOLATING EFFECTIVE LAGRANGIAN 

In this section, we lay out a general approach to calculating the neutron EDM using the 
effective Lagrangian method, independent of new physics. In this approach, one integrates 
out all heavy particles including SM gauge bosons and heavy-quarks. The resulting flavor 
neutral CP-violating effective Lagrangian has an expansion in terms of operators consisting 
of light-quark fields, u, d, and s and the gluon field G^ v ', with increasing dimensions, 

£C p-odd = J2 £f-°dd j (2) 

i=A,... 

where Li contains i-dimensional operator. The Wilson coefficient of each operator depends 
on the renormalization scale \i which in the end will be chosen as hadronic physics scale, 
about 1 GeV or lattice cut-off 1/a, where a is lattice spacing. The QCD operators also 
depend on the renormalization scale, but physical results do not. 
At the lowest dimension, there are two CP-odd operators, 

^ 4 = ~^r^ G ^ ( ^V + y^"W7 5 g , (3) 

OZ7T z * — ' 
Q 

Through SU(3) chiral rotations, the CP-odd quark-mass term can be rotated into a chiral 
singlet. Furthermore, one can eliminate either GG or the singlet quark-mass term through 
{7a (1) chiral rotation q — > e* 75 "g. 
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At dimension-five level, there are two kinds of flavor neutral P-odd and CP-odd operators, 
namely, the quark EDM operators and the chromo electric dipole moment (CDM) operators, 



A = E dE M° E M + E dC M°°M > ( 4 ) 

where = —^qcr^ijsqF^ and O^ 7 = — ^qcr^ v i^^t a qG a ^ u1 and and are the electro- 
magnetic and gluon field strengths, respectively, and t a are generators of the SU(3) gauge 
group. The one- loop evolution equations are [19J 

where 6/ = 11 — 2nj/3, rif is the number of quark flavors. It is easy to see that the 
dependence of the evolution of the quark CDM on nj is the same as that of the strong 
coupling, since they are both derived from wave function renormalization of the gluon field. 
At dimension- six, there are a number of four-quark flavor-neutral CP-odd operators and 



the Weinberg's three-gluon operator [23 



£ 6 = E^(A0°4*(/i) + C g (ji)O g (ji) , (7) 

i 

where the four-quark CP-odd operators can be divided into two groups. The first group 



includes operators with two different light flavors [24 



On = (gi75?)(?V) > 

0i2 = (qq)(qii5q') , 

On = (Vlst a q)(q't a q r ) , 
O22 = (qt a q)(q'ilst a q r ) , 

3 = (qii^ 1 'q)(q 'a , 

04 = (^^t a q)(q'a, u t a q') , (8) 
where q, q' — u,d,s and q 7^ q' '. The second group includes operators with one quark flavor 

o[ = {qiim){qq) , 

0' 2 = (qil,t a q)(qt a q) ■ (9) 



The Weinberg operator is defined as 



O g = - l -f^ e ^Gl p G b SG c ap , (10) 



where e 0123 = I. 

It is not difficult to see that all the dimension-six operators listed above are CP-odd. In the 
first group there are two different flavors in each operator. These operators are constructed 
by a pseudoscalar current coupled to a scalar one, a pseudo-tensor current coupled to a tensor 
one. No operator is constructed from an axial-vector current coupled to a vector current 
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since CP-odd operators cannot be generated in this way. And since 75 cr^ = f e^^cTo^, 
the two operators O3 and O4 are enough to describe the CP-odd pseudo-tensor and tensor 
coupling. Therefore, the first group includes all the P and CP-odd four-quark operators 
constructed by two different quark flavors. For the one-flavor case, the first four operators 
in the first group automatically become 0[ and 0' 2 in the second group, and using the 
Fierz transformation one can easily see that the operators described the tensor-pseudotensor 
coupling are not independent of 0[ and 0' 2 . Therefore, the second group includes all the 
flavor-neutral CP-odd four-quark operators with single quark flavor. 

The leading-order QCD evolution equations for dimension-six operator are as follows 
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The anomalous dimension of the Weinberg, 7 g9 , has been calculated in the literature [25 
jgg = — Ca/2 — Uf, where Ca = 3. The dimension-six operators mix with the dimension- 
five operators when scale evolves, however at the energy scale where only the light quarks 
exist, the mixing can be neglected because the dimension-five quark EDM and CDM are 
chirality flipping and thus proportional to the quark mass. At higher energies, the mixing 
is important and we will discuss it in the following sections. 

There is no mixing between the Weinberg operator and the four-quark operators listed 
in Eqs. (jHJ) and Q). To see this, we can decompose the four-quark operators into irreducible 
representations of the SU(3)l x SU( 3) n chiral group and only (3,3), (6,6), (8,8) and their 
conjugate representations are found [22|. On the other hand, the three-gluon operator is a 
chiral singlet. QCD evolution maintains the chiral structure of operators. 

When scale changes, the pure quark-gluon CP-odd operators generate perturbative con- 
tributions to quark EDM through the following T-product 



d A x T eA, 



X 



(14) 



where the summation neglects the quark EDM operator itself. The contributions are diver- 
gent so they induce additional running of the CP-odd operators. The contributions from the 
dimension-six operators are proportional to the mass of light quarks and can be neglected. 
The only large contribution is from the quark CDM operator, whose running has an effective 
inhomogenous term, 
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where Q q is the electric charge of the quarks and gs is the coupling of strong interaction. In- 
versely, the quark EDM operators can also generate quark CDM operators through the elec- 
tromagnetic interaction which is, however, proportional to the electromagnetic fine-structure 
constant. 

Therefore, omitting the ^-contribution, one can define the following electric dipole form 
factor 

= (iV(4)|J>>)Of(0;M) 



q 

*4„ 



+i / (fx T 



eA,(x)j? m (x) ^cf(/i)O g c (0;/i) 



\7(q)N(h)) , (16) 



where = k% — k± and U n is the wavefunction of neutron and e u is the polarization of 
the incoming photon. The static nEDM is just the zero-momentum limit of the form factor 

III. WILSON COEFFICIENTS IN LRSM 

Following the previous section, we make calculation of nEDM in the mLRSM by first 
evaluating the Wilson coefficients of the effective quark-gluon operators at the electroweak 
scale, and subsequently running them to hadronic scale. The detail of the model can be 
found in Ref. [15|, in which the spontaneous CP- violation is controlled by a phase angle a in 
the Higgs sector, and additional parameters of the model include, among others, the masses 
of the right-handed gauge boson and the new Higgs bosons. In the following subsections, we 
study the Wilson coefficients of various CP-violating operators separately. We will ignore 
the contribution of the #-term as it will usually generate a much too large nEDM: We assume 



certain mechanisms such as Peccei-Quinn symmetry [18( is in operation to suppress it. 



A. CP-Odd Four-Quark Operators 

To leading order, diagrams in Fig. [TJ generate the CP-odd four-quark operators induced 
by the exchange of gauge bosons and Higgs bosons. The operators are listed in Eq. (|Sj) and 
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(J9]). The corresponding Wilson coefficients can be easily read through the diagrams, 



C ab = ^^^2C\m{e- ia Vl b V^ b *) + ^^\m{C aa D bh " 



H2 



Cg = -^^ s i n 2(lm(e~ ia V^V^ b *) + ^^lm(C aa D bb " 

Ho 



H2 

C$ = V8G F sm2CIm(e- ia V£ b V* b *) 

1 V&Gf , 2 2\/- t / -ia-i ra0T ral3*\ 

+^72— ( m a -m b )£hn(e V L P V R P ) , 

IV1 H 2 

C% = -V8G F sin2(Im{e- ia V£ b V£ b *) 

- — (ml- ml)i Imie^Vfvf*) , (17) 




(b) 

FIG. 1: Effective four-quark operators generated by integrating out Wi-boson: (a) the diagrams 
in the full theory and (b) the effective operator. 
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Cf = ^(rnl-mlKMe^Vfvr) 

1V1 H 2 



2V8G 



C«f = ~ v M ~ F lm{C aa C a ' a '*) , 



l H 



2VSG 



C i2 = -^f^lm(C aa C a ' a '*) , 



Cg = -^^Im(IW' a '*) , (18) 

#0 



where a, a' e u,c,t, a ^ a' and b,b' E d,s,b, b 7^ b', C = V L M D V^ - 2£e ia M Ul D = 
VlMuVR—2^e~ ia M D , M Ho is the mass of the flavor changing neutral Higgs (FCNH) and M H2 



is the mass of which is a charged Higgs in mLRSM |15[]. Mu and Mo are diagonalized 
quark mass matrices. ( is the mixing angle between the lefthanded and righthanded W- 
bosons that 

. 4m b [MA 2 

sm2C~-r — , 19 

m t \M 2 J 

where r = (mt./mk)£ and ^ is the ratio between the two vevs of the Higgs bidoublet in 



mLRSM [15(. The contributions due to the Higgs exchanges are always proportional to 
quark masses. Since we are only interested in operators with at least two of the quarks 
being light, the Wilson coefficients are always proportional to at least one light quark mass, 
or they are proportional to heavy quark masses but must be suppressed by the non-diagonal 
CKM matrix elements. Furthermore, the mass of FCNH is strongly constrained to very large 
value by the mass differences and the CP-violating decay properties of the neutral K-bosons 



and B-bosons [la, [26| , and detailed calculation shows if 2 i s as heavy as FCNH. If we are 
interested in the case of a few TeV right-handed W-boson mass, we can safely neglect the 
Higgs exchange contributions. Then at the electroweak scale the Wilson coefficients of the 
CP-odd four-quark operators can be simplified to 

\/8C 

C{1 = -Cf 2 = - ^sin2C Im(e- ia V£ b VS b *) , 

6 

C£ = -C 2 t = v / 8G F sin2C Im(e- ta V£ b Vg b *) . (20) 
We will take this simple limit in the following discussion. 



B. Quark EDM and CDM Operators 

The one-loop contributions to the quark EDM from the gauge interactions are shown in 
Fig. [21 where the internal wavy lines represent the light charged gauge-boson W\ which is 
dominated by Wl, but has a small admixture of Wr. The dashed lines represent the charged- 
Goldstone boson present in Feynman gauge, and the external wavy line is the static electric 



S 



field or photon. Diagrams a) and b) have the photon interacting with the quarks directly, 
and these from c) to f) have the photon interacting with charged bosons. For the quark 
CDM case we have the first two diagrams only with the external wavy line representing a 
gluon. 

These diagrams have been calculated in the literature long ago our result is somewhat 
different from theirs in the infrared part. The CP-odd part of the diagrams in Fig. [2] can be 








FIG. 2: One-loop contribution to quark EDM. The internal wavy lines represent the W-boson 
contribution and the dashed lines the corresponding Goldstone bosons. 



expressed in terms of the coefficients of the EDM and CDM operators. For the up quark, 
we have (if Of + d^O^ with, 



d u = T^i Y\ m di eV^G F sin2Q\m{e- ia VtV^) 

107T Z z — ' 

i=d,s,b 



1 



X 



1 , 1 



4n + 3rf rf H — n In n rf kr; , . 

1 - n 3 V 3 1 3 J 2 2 1 1 ' ' 



d u = T^Yl rn A g a y/%G F ^20m.{e- ia V^Vt) 

i=d,s,b 

1 3 is 3 , 
x — 1 r,- r - H — r,Tnr,- 

1 - r,- 3 V 4 4 4 1 2 1 1 



(21) 
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And for the down quark, the contribution is d^Of + d^O^ with 
d d = J2 m ul eV8G F sm2Clm(e- ia Vl l V^*) 

i=u,c,t 

x (i_ r ()3 " ~J r * < ~ 12 ri ^ * ~ 2 ri 7 ' 
d C d = --^ V m m ( 7s v / 8G F sin2CIm(e- ia \/2 1 \4 1 *) 

l07T z / — ' 

x (T^)a( 1 -^-k + ^ lnr ' (22) 

It is easy to see that this mixing angle is suppressed by the mass ratio of bottom and top 
quarks and by the ratio of the left and right handed IF-boson masses. m U i are the masses 
of up-type intermediate quarks, = m di /M^, r[ = m^jMyy, Vl and Vr are the left and 
right-handed CKM mixing matrices, respetively, a is the spontaneous CP phase mentioned 
earlier. 




FIG. 3: Higgs-induced quark EDM. The dashed lines here represents the Higgs bosons. 



In mLRSM, H2 also gives contribution to the quark EDM and CDM. The relevant 
diagrams are shown in Fig. [3J and the result is 



dt 



]C 



E 

a£{d,s,b} 



167T 



2m a (m 2 u - ml) 



e Im(e-* a Vl a V^ ) 



3 — Arj + r 2 + 2 In Tj 
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9s V8G F 2 - ma{ml 



- 1 + r 2 — 2r j In rj 

2(-l +r J ) 3 
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Im(e-^^ L la ^ a *) / A, 



(23) 



d E d 



a d 



E_J_ A /o r r 2m a (m 2 m|) -ia T/ al T/ al*N 
16^ V8Gf M| * Im(e ^ ^ } 

3 - 4r' + rf + 2 In r' - 1 + r' 2 - 2r' In r' 
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a£{«,e,t} 



2(-l + r;.) 3 
9s V8G F 2 - ma{ml 



16tt 



^ 2 



2(-l + r;.) 3 

— ^(e-^v^ys 1 *)- 



4r'+rf + 21nr' , s 
J , (24) 



■1 +r]f 
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in which 



T ' ~ Mk • T > ~ M%, ' (25) 

Therefore, if the right-handed PL-boson has a moderate mass, say, a few TeV, the contri- 
bution from H2 to the quark EDM and CDM can be neglected in comparison to that from 
the right-handed gauge boson. 

Actually, there are both long-distance and short-distance contributions from the one-loop 
diagrams in Fig. [2] and Fig. [31 The short-distance contributions come from the integration 
region where the internal momentum is around M^; and the long-distance one from the 
loop momentum around the internal light quark masses. Due to asymptotic freedom of the 
strong interaction, the short- distance contributions can be calculated accurately using per- 
turbation theory. The long-distance contributions, however, suffer from non-perturbative 
QCD effects, and the only known way to calculate it correctly is by Lattice QCD. In the 
matching calculation, the long distance contribution has to be subtracted to obtain the Wil- 
son coefficients, which is shown in Fig. |H This contribution can be calculated using a certain 



FIG. 4: Long-distance contributions to quark EDM and CDM through CP-odd four-quark opera- 
tors. 



UV regulator, such as dimensional regulation or momentum cut-off or lattice regularization. 
Any regularization preserving a certain Fierz identity will give a zero answer as the loop in- 
tegral involves only the photon or gluon external momentum. Other regularizations, such as 
naive dimensional regularization, will find a finite contribution. One must be careful though 
that the vanishing of long-distance contribution is only true at one-loop level: as soon as 
the QCD corrections are taken into account, the result becomes non-zero. Therefore, to 
the leading order, we can directly read off the Wilson coefficients of quark EDM and CDM 
operators from Eqs. pT]) . (1221). and (1231. 



C. Weinberg Operator 

In mLRSM, the Weinberg operator can be induced from diagrams in Fig. Since the 
result is proportional to the quark masses, the leading contribution comes from the third 
generation of the quarks running in the loop. These are two-loop diagrams, the Weinberg 
operator comes out after one integrates out the internal quarks and bosons entirely. If one 
follows the effective theory approach, in which the top quark and the W-boson are first 
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FIG. 5: Diagrams contributing to Weinberg operator in mLRSM. The first diagram is induced 
by the W-boson exchange, the second by Goldstone exchange and the third by the charged Higgs 
boson. 



FIG. 6: Contribution to the three-gluon vertex after integrating out the top quark, the Higgs boson 
and the W-bosons. The black dot labels the bottom quark CDM operator. 



integrated out, the CDM operator of the bottom quark emerges and one can get its wilson 
coefficient from Eq. (|22|) . 

Then from Fig. [6j one gets the major contribution to the Weinberg three gluon vertex. 
Because this diagram would diverge quadratically in the infrared if the mass of the bottom 
quark was zero, this diagram should be proportional to However, chirality flipping 

is needed or otherwise the fermion loop will vanish, so the numerator of the diagram must 
be proportional to m b . Combining the two effects together, this diagram is proportional 
to /m b , where d^ is the bottom quark CDM which is proportional to m t . Therefore 
this diag ram has an enhancement of a factor of m t /mb, about 40, which was first found in 
Ref. [27J. Detailed calculation gives the Wilson coefficient 




This contribution is seemingly large, however, it is suppressed by a numerical factor, 1/(1 — 
mf/Mf) 3 ~ —0.02 in Eq. fl22|) . Therefore, the effect of the enhancement is totally canceled. 
Furthermore, the evolution also makes the contribution of this operator to be smaller at 
the low energy region [25]. Therefore, we safely neglect its contribution to nEDM in the 
following calculations. 





C g (m b ) 



g 2 s {m b )d%(m b ) 
167r 2 m b 



(26) 
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D. Wilson Coefficients at Hadronic Scale Through Leading-Order QCD Evolution 



The coefficient functions above, and hence the quark-gluon operators, are calculated at 
the high-energy electroweak scale, which is not yet useful for practical calculations. We are 
going to remedy this by running down the scale in the composite operator by including the 
leading logarithmic pQCD corrections. When we change the scale, dimension-six operators 
will mix with each other and generate dimension-five operators, and dimension-five operators 
will also mix with each other. The Wilson coefficients for CP-odd four-quark operators are 
shown to the leading order approximation in Eq. (1201) . From Eq. (iTTj) the renormalization 
group equations (RGE) keep this relation, and other CP-odd four-quark operators are not 
generated by the running. Then one can redefine the operators 

Oab r^ a b r} a b 

Of = Og - O* , (27) 

with the Wilson coefficients Of* = Cf\ and Cf 6 = Cff, respectively. Therefore, the RGEs of 
the Wilson coefficients of the dimension-six operators can be written as 

^Cf(^) = ^QV) , (28) 

which shows that C\ grows as the scale goes down, whereas C 2 does the opposite. 

The RGE of the quark CDM operators are a little bit complicated. For d quark and 
s quark CDM operators, as we discussed before, the c quark internal line gives a large 



contribution. Therefore, the RGEs of d and s quark CDM operators can be written as [19 



16vr 2 



{<ya + b f /2-5)d°(ri . (29) 



The Wilson coefficient of the up quark CDM operator is one order of magnitude smaller 
than that of the of d quark due to that m s /m c ~ 1/10. In the above formula, 731 = 5/2, 
732 = —1, 733 = —14/3, and 5 = —4 is the anomalous dimension of the quark mass. Detailed 
calculation gives, at m c , the relevant Wilson coefficients are 

Cf d ' s \m c ) = 3.0C^ d ' s) (M L ) , 
C^ d ' s \m c ) = 0.87C 2 {d ' s) (M L ) , 

iK) = 1.7^Cf d ' s \M L ) + 0.34^cf^(M L ) + 1.6dl(M L ) . (30) 



were is the mass of the SM W-boson. 

The CP-odd operators generate additional running of the quark EDM operators through 
the electromagnetic interaction. The RGE of the down quark EDM operator can be written 



as 



191 



d^ d (/ " } = "3 (16tt2)2 74lC 2 (M) - ~[^2-^ d d (M) - l^r>44 - S)d d (m) , (31) 
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where 741 = 16/3, 743 = 16/9, 744 = —16/3, and similarly for the strange quark. The RGE 
of the electromagnetic coupling e does not depend on the strong coupling constant g s up to 
one-loop, therefore, can be treated as a constant. At the charm quark mass scale, one can 
get 

CK) = ^(0mC« d ' s \M L )+0.34Cf> s \M L )) 

+0.17< S (M L )+0.83< S (M L ) . (32) 

which shows the explicit contributions from the running of the four-quark operators as well 
as CDM operators. 



IV. NEDM IN MLRSM AND CONSTRAINT ON LEFT-RIGHT SYMMETRY 
SCALE 

In this section, we carry out the last step of the nEDM calculation in mLRSM by in- 
corporating the neutron matrix elements of hadronic operators. We collect the state-of-art 
results in the literature and use them to constrain the parameters in mLRSM. We find that 
in order to satisfy the current experimental bound on nEDM and the data on kaon-decay 
parameter e, the right-handed gauge boson Wr might be as heavy as 10 ±3 TeV. This bound 
is far higher than the bound obtained previously from the kaon mass difference, making it 
difficult to discover left-right symmetry at LHC. 



A. Hadronic Matrix Elements 

The most difficult part in calculating nEDM is to estimate the hadronic matrix elements. 
In the literature, many different approaches, such as the SU(6) quark model, bag models, 
QCD sum rules, and chiral perturbation theory have been used to make estimations. In this 
subsection, we summarize the results and get some idea about their uncertainties. 



1. Contribution from Quark EDM 

In the SU(6) constituent quark model, the matrix elements of the quark tensor operators 
are simple and scale-independent 0, [1] , leading to 

d ( $ = -\d E u +\d E d . (33) 

Although it has been suggested that one should use the constituent quark masses in the 
formulas of quark EDM [8], this is incorrect from the point of view of factorization. 
In the parton quark model discussed in it was found, 

tfW = -0.508df + 0.746df - 0.226df . (34) 



From the QCD sum rules, one gets [28 



4? = (1 ± 0.5) x 0.7(-0.25df + df ) . (35) 
Different approximations are largely consistent. 
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2. Contribution from Quark CDM 

The contribution to nEDM from the quark CDM in the constituent quark model is 0] 

d -=tv d " + lv dC ^ ^ 

y 9s y 9 s 

where g s is the coupling of strong interaction at the energy scale where the model is ap- 
plicable. In this calculation, the authors assumed first that the neutron is composed of 
constituent quarks, and then treated the gluon field inside the neutron as a background, 
neglecting its kinetic energy. Therefore, Eq. (I36I) can only be seen as an order-of-magnitude 
estimate. 



Weinberg's naive dimensional analysis has also been used to estimate this contribution 11 

EM, 

d N~^ (0{l)d c u +0{l)d c d ) . (37) 



In Ref. 13], the authors used the chiral perturbation theory to calculate the singular part 



of the long distance contribution, 



d N ^—{d c u + d c d ) . (38) 
9s 



And finally, QCD sum rules analysis in Ref. |28|] gives 



d% ] = (1 ±0.5) x 2^f(0.5d£ + d£) , (39) 
9s 



where g s is the strong coupling constant at 1 GeV, about 2.5. 



3. Contribution from Weinberg Operator 

The contribution from the Weinberg's operator Oy/ can be estimated by Weinberg's naive 
dimensional analysis [23[, which is an order-of-magnitude estimate 

d$ ~ eMC g (/i)/47r w 100 MeV e C 9 (lGeV) , (40) 
where M = AirF^ ~ 1190 MeV and fi is the hadronic scale taking as 1 GeV. 



On the other hand, the estimate based on QCD sum rules gives [29| 

d$ ~ (10 - 30)MeV e C g (l GeV) , (41) 

which is considerably smaller. In any case, because of the small coefficient function, the 
Weinberg operator contribution can essentially be neglected. 

4- Contribution from Four-Quark Operators 

The hadronic matrix elements of the four-quark operators have been studied and reviewed 
in Ref. 22|. In this work we will take the results from that paper. 
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B. Numerical Results 



As discussed in Ref . [15| , combining with the kaon indirect CP- violation e parameter, one 
can use nEDM to get the most stringent lower bound on the mass of the right-handed W 
boson in the context of the mLRSM. In Ref. (l5| . the authors used naive factorization [l(| 
to estimate the contribution of four-quark operators. However, this method for baryons 
may not be valid even in the large- A^^ limit, and the uncertainty is unknown. Therefore, 
we have assumed a very large error on their matrix elements and the resulting constraint 
on the left-right symmetry scale is not very strong. In a dedicated study of these matrix 



elements [22] , we have gotten a much better understanding on their contribution. In Ref. [22 
the contribution of four-quark operators to nEDM was separated into two parts, the direct 
contribution and the meson-condensate contribution. For the direct contribution, quark 
models were employed to calculate the hadronic matrix elements, which is only an order-of- 
magnitude estimate be. However, for the meson-condensate contribution, the factorization 
method was used to calculate the meson matrix elements, which can be justified in the 
large-Ac limit. Since the meson-condensate contribution dominates over the direct one, 
we believe that we reached a factor-of-two accuracy in the matrix elements of four-quark 
operators. 
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0.00 0.05 0.10 0.15 0.20 



FIG. 7: nEDM contributed from operators, uij^udd (short dashed red line), ui^uss (long dashed 
green line), down quark EDM and CDM operators (solid blue line). 



In mLRSM, after neglecting the contributions from FCNH and the charged higgs boson 
exchange, nEDM depends only on three parameters, r, a, and Mw R , where a is the new 
source of CP- violation. Therefore, if a = 0, nEDM predicted by the mLRSM will be the 
same as that predicted by SM, about five orders of magnitude smaller than the upper bound 



given by the current experiment [30] . Whereas for e, there are two new contributions in 
mLRSM [l5|, the Dirac phase in the righthanded CKM matrix inherited from the lefthanded 
CKM matrix, and the spontaneous phase a. The new contribution from the Dirac phase 
is enhanced compared to the similar contribution in SM due to the chiral enhancement in 
the hadronic matrix element (see Ref. [2l[ for a good review). The contribution of the 
spontaneous CP-phase a must be adjusted to cancel the contribution of the Dirac phase. 
Therefore, in mLRSM there is a tension between nEDM and e that one cannot only adjust 
a to suppress all the new CP- violation sources, and a large Mw R is needed. As a result, 
nEDM and e together give a lower bound on Mw R - 
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FIG. 8: Constraints on the mass of Wr and the spontaneous CP-violating parameter a from the 
kaon decay parameter e (Mh = oo, red dots; Mh = 50 TeV, blue dots) and nEDM (green dots). 
For nEDM, we use the current experimental upper bound as the constraint and for e we use the 
criteria that the beyond-SM-physics contribution should not exceed 1/4 of the experimental value. 



In this new study, we use the QCD sum rules to estimate the contribution of the quark 



EDM and CDM operators, and use the results in Ref. |22| for the contribution of the four- 
quark operators. Fig. [7] shows the contributions to nEDM from different operators at fixed 
M\y R and r. The result from the Weinberg operator is too small to be included in the figure. 
It is clear that the contributions from four-quark operators are much larger than from quark 
EDM and CDM operators. One way to understand this is that in mLRSM the quark EDM 
and CDM operators are generated in the same way as the four-quark operators. The quark 
EDM and CDM operators are generated through diagrams in Fig. [2] and the four-quark 
operators are generated through diagrams in Fig. [IJ The Wilson coefficients roughly have 
the following relations 

df - ^4c 4 ■ £ * ^#C 4 , (42) 



where A and A' are two proportionality coefficients, C4 is the Wilson coefficient of certain 
four-quark operators. Take the down quark EDM as an example, A can be written as 
sin 2 9cm c /m u ~ 15, where 9c is the Cabibbo angle. From QCD sum rules, nEDM con- 
tributed by the down-quark EDM operator is approximately the down-quark EDM itself, 
whereas the nEDM contributed directly from the four-quark operator can be written as [22| 



i07T z 

where Bq ~ 2.2 GeV is related to SSB of the chiral symmetry. Since Bq 3> Am^, nEDM 
directly from the four-quark operator ui^ 5 udd is much larger than the contribution from the 
down quark EDM operator. Indeed, this is a common phenomenon in left-right models and 
two-Higgs-doublet models, where the quark EDM and CDM operators are always generated 
by the triangle diagrams in Fig. [2J and the internal lines are always quarks. In other types 
of new physics models, the internal lines can be other kind of fermions. For example, in 
super symmetric models, they can be gauginos, and in extra dimension models, they can 
be KK-fermions, where the above relation between quark EDM operators and four-quark 
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operators is no longer hold. In these models, quark EDM and CDM operators might be 
more important that four-quark operators. 



Using the matrix elements in Ref. 22], we calculate the constraint from the nEDM and 
kaon-decay parameter e on the allowed parameter space of mLRSM. The result is shown in 
Fig. [HI The allowed parameter region by the experimental upper bound on nEDM is shown 
as green dots. The constraints from e-parameter depends strongly on the mass of the FCNH 
in the theory. We have shown two possible values of M# , 50 TeV and oo for simplicity. We 
assume for e the new contribution should not exceed 1/4 of the experimental value. From 
Fig. [8] one can see that the lower bound for the Mw R from nEDM and e is around 10 TeV. 
If we assume a factor of 2 uncertainty on the hadronic matrix elements, the actual bound is 
10 ± 3 TeV. This will make a direct detection of the right-handed gauge boson very difficult 
at LHC if it exits. 



V. CONCLUSION 

In this paper, we have studied nEDM in mLRSM systematically by using effective field 
theory approach. The formula for calculating nEDM is given in Eq. ffl6l) . The contribution 
of four-quark operators is found to be the most important. The contribution of Weinberg 
operator to nEDM has been discussed systematically. A numerical suppression is found 
which counteracts the infrared enhancement and makes the contribution of this operator 
negligible. We have found a lower bound on the mass of Wr which is about (10 ± 3) TeV. 



This is higher than what have been found before [15|, |3jj and certainly cannot be detected 
at LHC. 

In a more complicated non-supersymmetric scenario of LRSM, although the CP-violation 
pattern in the Higgs sector might be change, the tension between e and nEDM discussed 
in Sec. IV still exists. Therefore, one can also use this analysis to set a lower bound on 
the righthanded scale. In the supersymmetric LRSM, there are new CP-violation sources 
from the soft terms, which can contribute to both nEDM and e. Furthermore, in supersym- 



metric LRSM [32j, the lefthanded and righthanded CKM matrices must be equal to each 



other up to a sign, therefore, if one assumes certain scenarios of the breaking mechanism of 



supersymmetry, e itself can give a constraint on the righthanded scale [33 
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